Vortex Lattice Structure of Fulde-Ferrell-Larkin-Ovchinnikov Superconductors 
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In superconductors with singlet pairing, the inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state is expected to be stabilized by a large Zeeman splitting. We develop an efficient 
method to evaluate the Landau- Ginzburg free energies of FFLO-state vortex lattices and use it to 
simplify the considerations that determine the optimal vortex configuration at different points in the 
phase diagram. We demonstrate that the order parameter spatial profile is completely determined, 
up to a uniform translation, by its Landau level index n and the vortex Lattice structure and derive 
an explicit expression for the order parameter spatial profile that can be used to determine n from 
experimental data. 



I. INTRODUCTION 

Forty years ago Fulde and Ferrelli, and Larkin and 
Ovchinnikov 2 proposed that a superconducting state 
with an inhomogeneous order parameter would become 
stable when a singlet superconductor is disturbed by Zee- 
man splitting between electrons with opposite spin ori- 
entations. The Zeeman splitting can be due to either 
an external magnetic field or an internal exchange field. 
In spite of ongoing theoretical interest in the Fuldc- 
Ferrell-Larkin-Ovchinnikov(FFLO) state, evidence for its 
occurrence has been sparse until recently. The situation 
has been changed by experimental results suggestive of 
the FFLO state in organic, heavy-fermion, and possi- 
bly other unconventional supercondcutors^*Si2i2iiSiii. 
Recent experimental results on the heavy fermion com- 
pound CeCoIns, a quasi- two-dimensional d-wave super- 
conductor, are particularly compelling 910-11 . Indepen- 
dently, speculations that the FFLO state may also be 
realized in high density quark matter have caught the 
attention of the particle physics community, adding to 
the flurry of interest in this uncommon superconducting 
state. 

To establish the presence of the FFLO state unambigu- 
ously one needs to perform a phase sensitive experiment 
that detects order parameter oscillations and measures 
their wavelength. One possible phase sensitive probe 
based on the Josephson effect between a FFLO supercon- 
ductor and a BCS superconductor— was proposed some 
time ago by Yang and Agterberg. Although conceptually 
straightforward this experiment is yet to be carried out, 
possibly due to technical difficulties in its implementa- 
tion. More recently, there has been interest in the vortex 
lattice structure (VLS) of an FFLO state subjected to an 
orbital magnetic field, and particularly in the possibility 
of using the character of the VLS to infer the presence 
of an underlying FFLO superconductor 13 i 14 i 15 i 16 i 17 . The 
key ideai& is that, depending on the interplay between 
the orbital and Zeeman effects of the magnetic field, the 
order parameter of a FFLO state near its upper critical 



field can correspond to a high Landau level (LL) index 
Cooper pair wave function. Recent worlsiiiiiiSiiSiii on 
the FFLO VLS in specific situations has demonstrated 
that these high LL index VLS's can be very different 
from the triangular lattice Abrikosov VLS favored by low- 
est Landau level (LL) Cooper pairs. The case of FFLO 
states in an orbital field is the only example known to 
us in which bosons (Cooper pairs in the present case) are 
placed in high LLs; the absence of a fermionic Pauli exclu- 
sion principle implies that one cannot force bosons to oc- 
cupy a high LL simply by increasing the particle density. 
Electronic states in high LLs have recently been found 
to exhibit physics^ that can be quite different from that 
found in the lowest LL, suggesting that surprises might 
also occur in high LL bosonic physics. 

In this paper we apply techniques developed in studies 
of electronic states in a strong magnetic field (primar- 
ily motivated by the quantum Hall effect—) to study 
the vortex states of FFLO superconductors. Using 
these techniques, we develop a very efficient method 
to determine the VLS implied by a given Ginsburg- 
Landau free energy functional. As established by previ- 
ous studiesi^iii^iiSiii, the VLS is highly non-universal; 
we argue here that the VLS approaches a universal limit 
near the tricritical point, in the sense that the structure 
near this point in the phase diagram depends on the LL 
index of the order parameter only. We determine the 
VLS in this limit explicitly. We demonstrate that the 
order parameter spatial profile is completely determined, 
up to a uniform translation, by its Landau level index n 
and by the VLS, and derive an explicit expression for the 
order parameter spatial profile that can be used to deter- 
mine n from experimental data. If the LL index can be 
determined by comparing experimental and theoretical 
VLSs, the order parameter momentum in the absence of 
a field can also be determined. 

We lay out the theoretical method we have developed 
in the following section, and present the VLSs suggested 
by its application in section III. These techniques may be 
useful in addressing other aspects of high Landau level 
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bosonic physics in FFLO superconductors or elsewhere. 
In section IV we summarize our results and make a com- 
parison with related recent work. 

II. THEORETICAL METHOD 

To address FFLO VLS physics we follow a Ginzburg- 
Landau approach similar to that of Refs ll5llT^ . The free- 
energy density appropriate for the FFLO state, has the 
formal: 

F oc |(-V 2 - q 2 )^>\ 2 + a\i>\ 2 + 6|^| 4 + c#| 2 |W>| 2 
+ d[(V/) 2 (V^) 2 +^ 2 (V^) 2 ]+e^| 6 , (1) 

where a,b,c,d,e and q are parameters that depend on 
both temperature and Zeeman splitting. The free energy 
of Eq.Q is appropriate for an s-wave superconductor 
with completely isotropic gap and band dispersion, the 
case we will consider here. Anisotropy in either pair- 
ing interaction or single electron dispersion will lead to 
additional terms^. The fundamental difference between 
FFLO and BCS superconductors is expressed by the first 
term in F which describes the kinetic energy of the or- 
der parameter; in an FFLO superconductor this term is 
minimized when the order parameter carries a finite wave 
vector (or momentum) q; in a BCS superconductor q = 0. 
Thus far we have only taken into account the Zeeman ef- 
fect of the external magnetic field, which is appropriate 
for the case of a two-dimensional superconductor with 
an external magnetic field parallel to the plane. For a 
3D superconductor or for a 2D superconductor with the 
field B tilted out of system plane, orbital coupling must 
be accounted for by performing a minimal substitution 
W -»■ = (V - 2ieA/c)^ with V x A = B±z. It 
is the 2D case that we have in mind here because it al- 
lows the Zeeman field, which controls q and the other 
LG model parameters, and the field B± to be controlled 
independently. 

To find the order parameter ip that minimizes the free 
energy of the system near a continuous phase transition, 
one first determines the Cooper pair Landau level index 
that minimizes the coefficient of the quadratic terms in 
Eq. QJ. The eigenvalues of D 2 = (V — 2ieA/c) 2 are 
— (2n + l)/^ 2 , where I = \Jficj2eB is the Cooper pair 
magnetic length and n = 0, 1, 2, • • • is the LL index. For 
a BCS superconductor the kinetic energy is minimized by 
n = 0, i.e., ip is a lowest LL wave function. For a FFLO 
superconductor, however, the kinetic energy is minimized 
by the index n that minimizes \(2n + l)/£ 2 — q 2 \. The 
optimal value of n will be large when B± is small^ 

Because of the macroscopic degeneracy of a Landau 
level, the order parameter retains considerable freedom 
at Bj_ =/= even when n is fixed. As we explain below, this 
freedom is resolved by minimizing the remaining terms of 
Eq. Q with a properly normalized wave function ip(r); 
this corresponds to optimizing the vortex lattice. The 
overall magnitude of the order parameter is determined 



by balancing quadratic and higher order terms in the 
GL energy functional. For a BCS superconductor that 
is described by an ordinary Ginsburg-Landau theory for 
which n = and in which the |^| 4 term is the only addi- 
tion term, it is known that the optimal VLS is triangular. 
In the following we use techniques developed in studies 
of the quantum Hall effect to derive expressions for and 
understand the relationships between the higher order 
terms present in Eq. QJ. 

We assume that the VLS is that of a Bravais lattice so 
that there is one flux quantum per unit cell. The central 
quantity in our calculation is the Fourier transform of 
^(r)! 2 : / cPr|V'(r)| 2 e~' ik ' r . Since |VK r )| 2 is periodic, this 
integral is zero unless k = G, where G is a reciprocal 
lattice vector of the vortex lattice. For k = G, we define 

A G = / d 2 r|V(r)| 2 e - tG - r = (V|e- ,;G - r ^>- (2) 

J u.c. 

Here J c stands for integration within a unit cell of the 
vortex lattice, and we assume that ip(r) is normalized 
within a unit cell so that 

f d 2 i#(r)| 2 H^)=l. ( 3 ) 

J u.c. 

To proceed, we introduce the guiding center coordinate 23 

R = r + < 2 zx (-iV + (2e/c)A), (4) 

and separate r into r = R+(r— R). It is knownSSthat the 
guiding center R has no matrix element between states 
in different LLs, while (r — R) changes the LL index of 
a state without changing its guiding center dependence; 
as a consequence components of R and (r — R) commute 
with each othesSi. 

Two key observations simplify the calculations below. 
First \ip) is an eigenket of e _lGR . It follows that A G is 
the eigenvalue of this operator and, since the operator is 
unitary, must be a phase factor. As we discuss later the 
phase factor can be fixed by a convenient choice of spatial 
origin. The operator exp(— iG-R) results in a translation 
by G^ 2 x z in real space, which is a direct lattice vector, 
and multiplication by a phase factor that is gauge but 
not position dependent. The most straightforward way 
to prove this property is to construct ip(r) explicitly in a 
convenient gauge (say the Landau gauge), show that it 
is true in this gauge, and that an arbitrary gauge trans- 
formation (on ^i(r) and e _lG R simultaneously) does not 
change the result. The second key observations follows 
from the fact that r — R does not act on guiding center 
coordinates: 

Me- iG -< r - R >|V) =L n {G 2 f/2)e- GH2 / i (5) 

for any normalized n-th Landau level state. (Here L n (x) 
is the Laguerre polynomial.) 

It follows from the above properties that 

A G = (V|e- 4G '» = (V-|e- lG - R |^)^|e^ G - (r - R) |^) 
= A G L n (G 2 l 2 /2)e~ G2 ' 2 /\ (6) 
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with |Ag| = 1. Eq. © allows for a very straightforward 
calculation of the Abrikosov factor 



(n) 



A / d 2 r\ip{r) 

J U.C. 

J2lL n (G 2 e 2 /2)} 



G 

2 -G 2 ( 2 /2 



(7) 



G 



where A = 2ir£ 2 is the area of the unit cell which contains 
one flux quantum. In Eq. J7J the VLS is characterized by 
its reciprocal lattice vectors G. For BCS superconduc- 
tors described by ordinary Ginsburg-Landau free energy 
functional, the VLS is determined by minimizing /3 A °^ (24 
it is known that in the lowest LL (n = 0), is min- 
imized by the triangular lattice with = 1.1596j24 a 
result easily reproduced using Eq. Q, with the choice 

G = mgi + ng2, where gi = [\J 4ir/V3/£)x and g 2 = 

{\J ir/y/3/£)x + [\J \/3ir/£)y are the basis vectors of the 
(reciprocal) triangular lattice. For the square lattice, we 
have gi = (V2tt/£)x and g2 — {\/2tt / £)y , which yields 
(3^ = 1.1803, also in agreement with previously known 
results H 

To determine the VLS of a FFLO superconductor one 
needs to take into account the fact that the LL index n 
is not necessarily zero and the presence of contributions 
from terms other than |^(r)| in the free energy func- 
tional (JIJ. These can be done straightforwardly using 
extensions of Eq. Q that also follow from the sepa- 
ration of guiding-center dependence in the appropriate 
operators. For example the contribution from the \tp\ 
term is parametrized by 



(n) 

Ta 



A 2 f d 2 r|^(r)| 6 = Y, A Gl A G2 A_ Gl _ G 



G1G2 



= 2J A Gl A G2 A_ Gl _ G2 e" 



-(G 2 1 +Gl+G 1 -G 2 )e 2 /2 



G1G2 



x L n {G\£ 2 '/2)L n {G 2 2 £ 2 : /2)Z„(|Gi + G 2 \ 2 £ 2 /2) 

= Y cos ^ ■ ( G i x Q 2 )£ 2 /2]e-^ +G2+G ^ 12 ' 2 



G1G2 



x L n (G 2 1 £ 2 /2)L n (G 2 2 £ 2 /2)L n (\G 1 + G 2 \ 2 £ 2 /2), (8) 
where we used the property 



-i(Gi+G 2 )-R _ -1G1R -iG 2 'R !z (GixG 2 )£ 2 /2 



(9) 



which follows from the commutator between different 
components of R-i [R x , R y ] — i£ 2 . It follows that 

A_ Gl _ G2 = A-iA-y^xG,)^. (io) 

Similarly, using the fact that ip is an eigenfunction of 
the operator (V — 2ieA/c) 2 , the contributions from the c 
and d terms in Eq. JIJ can be expressed as combinations 
of (3^ ancU£ 



1 11 — 



d 2 r|^(r)| 2 |D^(r)| 5 



(11) 



To evaluate I 22 \ we introduce 

A G = / d 2 r|D^(r)| 2 e - 4G ' r 

d 2 n/>*D[D^(r)e- iG - r ] 

= - f d 2 r^*[D 2 V(r) - i(G • D)i/.(r)]e~ lG tl2) 

It is easy to show using the fact that tp is an eigen wave- 
function of D 2 that the second term in the final line of 
Eq. ijT^t does not contribute to the integral. We thus 
find 



A G = {2n+l)A G /£ 2 



(13) 



and 



1$ = Y A G A_ G = (2n + l)^/£ 2 . (14) 



We have thus found that I 22 ^ is proportional to (3^ in 
particular, they have exactly the same dependence on 
the VLS. This fact simplifies the determination of the 
optimal VLS significantly. 



III. RESULTS 

In tables I and II we have listed the VLS that mini- 
mizes (3^ and 7^ respectively, for the order parameter 
in the nth Landau level. Fig. contains plots for small 
Landau level indices of the square of the order parame- 
ter magnitude for the vortex lattices that minimize (3^ . 
For n = we find that that both /3a and 7^4 are mini- 
mized by the triangular lattice, in agreement with previ- 
ously known result o 15 ' . Evidently the situation changes 
drastically for n =/= 0. The competition between different 
potential VLS's can be understood in terms of the sim- 
ple reciprocal lattice expression of j3^\ Eq. Q. For 
n = 0, the Laguerre polynomial is L n —o(x) = 1, so it is 
clear that (3a is minimized by placing the reciprocal lat- 
tice vectors (G's) as far away from the origin as possible 
given the fixed unit cell area, which is what a triangu- 
lar lattice does. For n > 1, on the other hand, L n (x) is 
a polynomial with nodes whose positions depend on n; 
one may thus lower (3^ by, for example, choosing G's 
such that L n (G 2 £ 2 /2) is close to zero for some reciprocal 
lattice vectors, leading to very non-standard VLSs under 
most circumstances. Similar considerations apply to the 
optimal VLS that minimizes 7^- Overall, the optimal 

values of both j3^ and 7^ increase with n, reflecting 
the fact that the order parameter wave function becomes 
increasingly oscillatory (and therefore less uniform) as n 
increases. For n = \ and 3 VLS's are characterized by or- 
der parameter concentration in widely spaced stripe-like 
regions with rapid but weak longitudinal modulations. 
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FIG. 1: (Colo? online) Plots of the square of order parameter magnitude that minimizes p£ , for Landau levels < n < 6. 
The lines (when present) indicate the region of a unit cell of the vortex lattice. Lengths are in unit of the magnetic length I. 
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TABLE I: The optimal vortex lattice structure that minimizes 
P^T in different Landau levels. The vortex lattice structure 
is parametrized by the aspect ratio between the two basis 
vectors, and the angle between them. 
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TABLE II: The optimal vortex lattice structure that mini- 
mizes 7^ in different Landau levels. The vortex lattice struc- 
ture is parametrized by the aspect ratio between the two basis 
vectors, and the angle between them. 



The stripes are characterized by lines with definite guid- 
ing centers and odd parity perpendicular quantum wave 
functions, yielding a series of lines along which the order 
parameter is small. For even n and for larger n the op- 
timal VLS's tend to be more close-packed, with lattices 
that are close to square occurring frequently. 

Since both the quartic terms (6, c and d terms in Eq. 



Q) and the sixth order term (the e term) contribute to 
the total free energy, and f3^ and 7^ depend on the 
lattice structure in different ways for n > 1, the optimal 
VLS will depends on an interplay between the two, and 
is therefore non-universal; this is very different from the 
Abrikosov lattice of the BCS superconductor, which is al- 
ways triangular. In the following we discuss two limiting 
cases in which the optimal VLS approaches a universal 
limit, in the sense that it depends on the Landau level 
index n only. 

(i) It is knownSi that at the tricritical point which oc- 
curs for T ss 0.56T C (the maximum temperature at which 
the FFLO state is stable) both the order parameter wave 
vector q and the coefficient of the quartic term b vanish. 
Near this tricritical point, the contributions of the quartic 
terms to the free energy vanishes (the terms proportional 
to c and d both involve gradients are proportional to q 2 ). 
The optimal VLS is therefore determined by minimizing 

(n) 

7^ alone and is given by the entries listed in table II. 

(ii) Along the H C 2 line, the transition can be cither 
first order or second order. The transition is second 
order when the combined contribution to the free en- 
ergy from the quartic terms is positive, in which case 
the contribution from sixth order term is negligible near 
the phase boundary, due to the smallness of the order 
paramter. As shown above, the contribution to the free 

(n) 

energy from all the quartic terms is proportional to f3 A 
with lattice structure independent prefactors. Minimizing 
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the free energy is therefore equivalent to minimizing (3^ i 
resulting in an optimal VLS that depends on n only and 
given in this case by the entries listed in table I. There 
is, however, an important caveat in this case since the 
conclusion is reached based on the free energy Eq. 
in which quartic terms that involve more gradients are 
neglected; these terms are in principle comparable to the 
terms kept, because the order parameter wavevector q is 
finite, and it may not be possible to express their contri- 
butions in terms of (3^ alone. These terms do not affect 
case (i) however because q approaches zero there. 

Once the VLS is determined, one can obtain the de- 
tailed spatial dependence of the order parameter in the 
nth LL by counter- Fourier transform Eq. (J2J): 

= rJpE A L„(GV/2)e- G *'''V G "(15) 

G 

By parameterizing the reciprocal lattice vector G = 
TOigi + m 2 g 2 where gi and g 2 are the two basis vec- 
tors of the reciprocal lattice, while mj and m 2 are two 
integers, we find 

Ag = (V , |e~^ miSl+ " l2S2 - ) ' R |'i/') = e "" im2£ '(sixg2)^ 2 /2 
x (^|e _imigl ' R e _ " n2S2 ' R |^) = (— l) m l m 2e i (" ^leI+m^ 'i 3 j , 

where e lBl and e 102 are the eigen values of e~ Jgl ' R and 
e -jg 2 -R respectively; we may choose 0\ = 9 2 = by se- 
lecting the spatial origin in the unit cell. We thus obtain 

l^(r)| 2 = ^ £ (-l)™L n (|m lgl +m 2g2 | 2 £ 2 /2) 

77117712 

x e - |miSl+m2S2|2 ^ /4 cos[(migi + m 2 g 2 ) • r](17) 

Experimentally one can determine the VLS (and thus 
gi and g 2 ), by neutron or muon scattering for exam- 
ple, while the magnitude of the local order parameter (or 
gap) may be measured by STM. By comparing the mea- 
sured results with Eq. (|17|l one may be able to determine 
the LL index n. This allows for an estimate of the or- 
der parameter wavevector in the absence of the orbital 
magnetic field: q w \j2n + T/£, which is an alternative to 
measuring q using the Josephson effecti^. It is important 
however to keep the following caveats in mind when mak- 
ing this comparison, (i) We have assumed here that the 
order parameter wave function is in a given Landau level 
exclusively. This is true only when one is very close to the 
portions of i? c2 line where the transition is second order, 
or to the tricritical point mentioned above. Away from 
these regions mixing of different Landau levels becomes 
increasingly important, and this can affect our results 
significantly, (ii) Our analysis is based on a Ginsburg- 
Landau free energy that is completely isotropic. Real sys- 
tems often have some anisotropy, due to either the pair- 
ing interaction (as in a d-wave superconductor) or single- 
electron dispersion (band structure effect). It should be 



possible however to generalize our methods to anisotropic 
cases. 



IV. SUMMARY 

In this work we have presented and applied an effi- 
cient method to determine the vortex lattice structure 
(VLS) of an FFLO superconductor, based on an appro- 
priate Ginsburg-Landau free energy functional. In our 
approach all calculations are performed in a completely 
gauge independent manner. We have shown that the VLS 
is universal in certain limiting cases, in the sense that it 
depends on the Landau level index of the order parameter 
only. Independent of the VLS, which is not universal in 
the general case, our results enable an experimental de- 
termination of the Landau level index which follows from 
the constraints on spatial structure that exist when only 
guiding center degrees of freedom are available. Given 
the Landau level index, one can estimate the order pa- 
rameter wavevector in the absence of orbital magnetic 
field. 

Fundamentally, our approach is similar to that used 
in Refs. IT5IIT6I in that it is based on minimizing an ap- 
propriate Ginsburg-Landau free energy functional. The 
efficiency of the new calculational technique, however, 
allows one to determine the vortex lattice structure for 
a wide range of parameters, instead of only the specific 
cases studied earlier. Our methods can be applied to any 
problem in which bosons occupy high Landau levels. In 
addition to vortex states of FFLO superconductors, an- 
other possibility is bosonic cold atoms in a rotating trap 
with an optical lattice potential that has been engineered 
to have the minimum band energy away from the zone 
center. 

In a recent paper^ Klein studied the VLS of FFLO 
superconductors starting directly from a weak-coupling 
BCS model (which is a generalization of earlier work by 
Eilenberger— ) , instead of using a Ginsburg-Landau free 
energy functional. In his work the VLS was determined 
solely from the terms that are quartic in the order pa- 
rameter in the free energy, while as we see here terms 
beyond quartic can be crucial in certain regimes. On the 
other hand it appears that he was able to evaluate the 
quartic terms without using a gradient expansion, which 
appears to be an improvement on the usual Ginsburg- 
Landau approach; this might be important under some 
circumstances. We note that we find the dependence of 
VLS on the Landau level index n appears to be quite 
complicated, without any obvious indication of the sim- 
plification in the limit n ~ > 00, suggested in Ref. Il7l 

Acknowledgments 

We thank Q. Cui for technical assistance. KY was sup- 
ported by NSF grant No. DMR-0225698, and by the Re- 
search Corporation. AHM was supported by NSF grant 



6 



No. DMR-01 15947 and by the Welch Foundation. 



1 P. Fulde and A. Ferrell, Phys. Rev. 135, A550 (1964). 

2 A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 20, 
762 (1965). 

3 For a review, see R. Casalbuoni and G. Nardulli, Rev. Mod. 
Phys. 76, 263 (2004). 

4 K. Gloos, R. Modler, H. Schimanski, C. D. Bredl, C. 
Geibel, F. Steglich, A. I. Buzdin, N. Sato, and T. Ko- 
matsubara, Phys. Rev. Lett. 70, 501 (1993). 

5 R. Modler, P. Gegenwart, M. Lang, M. Deppe, M. Wei- 
den, T. Lhmann, C. Geibel, F. Steglich, C. Paulsen, J. L. 
Tholence, N. Sato, T. Komatsubara, Y. Onuki, M. Tachiki 
and S. Takahashi, Phys. Rev. Lett. 76, 1292 (1996). 

6 M. Tachiki, S. Takahashi, P. Gegebwart, M. Weiden, M. 
Lang, C. Geibel, F. Steglich, R. Modler, C. Paulsen, and 
Y. Onuki, Z. Phys. B 100, 369 (1996). 

7 P. Gegebwart, M. Deppe, M. Koppen, F. Kromer, M. Lang, 
R. Modler, M. Weiden, C. Geibel, F. Steglich, T. Fukase, 
and N. Toyota, Ann. Physik 5, 307 (1996). 

8 J. Singleton, J. A. Symington, M.-S. Nam, A. Ardavan, M. 
Kurmoo, and P. Days, J. Phys. Condens. Matter 12, L641 
(2000). 

9 H.A. Radovan, N.A. Fortune, T.P. Murphy, S.T. Han- 
nahs, E.C. Palm, SW. Tozer, and D. Hall, Nature 425, 
51 (2003). 

10 A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and 
J. L. Sarrao, Phys. Rev. Lett. 91, 187004 (2003). 

11 C. Martin, C. C. Agosta, S. W. Tozer, H. A. Rado van, E . 
C. Palm, T. P. Murphy, J. L. Sarrao, cond-mat /0309125 

12 K. Yang and D. F. Agterberg, Phys. Rev. Lett. 84, 4970 
(2000). A related experiment was proposed in L. Bu- 
laevskii, A. Buzdin and M. Maley, Phys. Rev. Lett. 90, 
067003 (2003). 

13 H. Shimahara and D. Rainer, J. Phys. Soc. Jpn. 66, 3591 
(1997). 



1 U. Klein, D. Rainer, and H. Shimahara, J. Low. Temp. 
Phys. 118, 91 (2000). 

15 M. Houzet and A. Buzdin, Europhys. Lett. 50, 375 (2000). 

16 M. Houzet, A. Buzdin, L. Bulaevskii, and M. Maley, Phys. 
Rev. Lett. 88, 227001 (2002). 

17 U. Klein, Phys. Rev. B 69, 134518 (2004). 

18 L. N. Bulaevskii, Sov. Phys. JETP 38, 634 (1974). 

19 M. M. Fogler, A. A. Koulakov, and B. I. Shklovskii, Phys. 
Rev. B 54, 1853 (1996); R. Moessner and J. T. Chalker, 
Phys. Rev. B 54, 5006 (1996); M. P. Lilly, K. B. Cooper, 
J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Phys. 
Rev. Lett. 82, 394 (1999). R. R. Du, D. C. Tsui, H. L. 
Stormer, L. N. Pleiffer, K. W. Baldwin, and K. W. West, 
Solid State Commun. 109, 389 (1999); E. H. Rezayi, F. 
D. M. Haldane, and K. Yang, Phys. Rev. Lett. 83, 1219 
(1999); F.D.M. Haldane, E H. Rezayi, and K. Yang, Phys. 
Rev. Lett. 85, 5396 (2000). 

20 The Quantum Hall Effect, 2nd Ed., edited by R. E. Prange 
and S. M. Girvin (Springer, New York, 1990). 

21 A. I. Buzdin and H. Kachkachi, Phys. Lett. A 225, 341 
(1997). 

22 D. F. Agterberg and K. Yang, J. Phys.: Condens. Matter 
13, 9259 (2001). 

23 See, e.g., F. D. M. Haldane, in Rel|U A.H. MacDonald in 
Proceedings of the Les Houches Summer School on Meso- 
scopic P%sicsElsevierAmsterdaml995edited by E. Akker- 
mans, G. Montambeaux, and J.-L. Pichard. 

24 See, e.g., D. Saint-James, G. Sarma and E. J. Thomas 
(eds.), Type II Superconductors, Pergamon Press, New 
York (1969). 

25 G. Eilenberger, Z. Physik 190, 142 (1966); Phys. Rev. B 
153, 584 (1967). 



